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FRAGMENTATION AND COLLAPSE OF TURBULENT MOLECULAR CLOUDS 
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Abstract 

We performed simulations of self-gravitating hydro- 
dynamic turbulence to model the formation of filaments, 
clumps and cores in molecular clouds. We find that when 
the mass on the initial computational grid is comparable 
to the Jeans mass, turbulent pressure is able to prevent 
gravitational collapse. When the turbulence has damped 
away sufficiently, gravitational collapse can occur, and the 
resulting structure closely resembles the pre-singularity 
collapse of an isothermal sphere of Pcnston (1969). If 



several Jeans masses are initially placed on the grid, tur- 
bulence may not be sufficient to prevent collapse before 
turbulence can be significantly damped. In this case, the 
cores have density structures which are considerably shal- 
lower than expected for an isothermal gas, and resemble 
the solutions for a logatropic equation of state. 
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1. Background 

There is abundant evidence that stars form in clumps 
within molecular clouds. In order to understand the mech- 
anism by which stars form, we therefore must understand 
the environment within a molecular cloud. Observations 
show that molecular clouds are composed of a hierarchi- 



cal network of filaments, clumps and cores ( Mizuno et al 
1995|; iJohnstonc fc Bally 1999| ; pahm et al. 2002| ). The 



most popular explanation for this substructure is that 



of turbulent fragmentation (Klessen & Burkert 2000; Os 



triker et al. 2001). In this scenario, turbulent fluctuations 
within the molecular cloud can compress the gas through 
shocks to high enough densities that the self-gravity of 
such regions becomes strong enough to induce collapse, 
resulting in a cloud which has been broken up into sev- 
eral condensing cores. The cores, due to their origins as 
fragmenting shocks, are arranged along filamentary struc- 
tures; the intersections of these filaments have many of 



these smaller cores attracted to their higher average den- 
sity, resulting in clusters which resemble the observed 
clumps. 

We performed simulations to test whether the struc- 
ture of cores produced in this turbulent fragmentation 
scenario are consistent with observation. 

2. Simulations 

Our simulations were performed using the NCSA ZEUS- 
MP hydrodynamics code. ZEUS-MP is a mesh-based Eu- 
lerian hydromagnetic code capable of performing 3D sim- 
ulations. The self-gravity of the fluid is calculated using 
a FFT gravity solver. Our simulations are stopped when 



the numerical gravitational instability criterion of Tru- 
love et al. (1997) is violated (Lj < 4L pix ). We use an 



isothermal equation of state and periodic boundary con- 
ditions. 

We start with a uniform initial density and no mag- 
netic field, at a resolution of 128 3 grid cells. Into this box 
we inject an initial series of plane waves, with their am- 
plitudes described by a Kolmogorov scaling law (Ef. cx 
fc~ 5 / 3 ) and with random phases and directions. We en- 
sure that V • v = 0. Our simulations can be described 
by two parameters, the RMS Mach number used to nor- 
malize the velocity fluctuations, and the initial number of 
Jeans masses on the grid, used to normalize the density. 
We chose our initial number of Jeans masses from the 



clump data of Lada et al. 1991 



3. Results 

We show the projected density of two of our simulations 
in Figure [T| each o f which is at the point where the Tru- 



love et al. (1997) condition is violated. In the simulation 



with an initial mass of 1.3Mj (where Mj is the initial 
Jeans mass), the relatively weak gravitational force was 
not able to induce collapse until much of the turbulent 
kinetic energy had dissipated. The resulting core is much 
more homogeneous than might be expected for a turbu- 
lent process. The density structure of this core is slightly 
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Figure 1. Simulation results at the time when the Truelove et al. condition is violated. Projected density is shown. 
Left: A run with initial Mrms = 5, nj = 1.3. The simulation has run for 2 crossing times (7.1 free-fall times). Right: 
A run with initial Mrms = 5.1, nj = 5. The simulation has run for 0.23 crossing times (1.28 free-fall times). 



Density Profile Evolution for a Clump Density Profiles of Many Clumps 




Figure 2. Density structure of cores. Left: A run with initial MrmS = 5, nj = 1.3. The density structure of the largest 
core is shown, from bottom to top, at 1.6, 1.7, 1.8, 1.9, and 2.0 crossing times. Right: Several cores for the run with 
initial Mrms = 5, nj = 5, after 0.23 crossing times. 



shallower than a Bonnor-Ebert sphere (Figure ||), sug- 
gesting that some small amount of turbulent pressure is 
providing additional support to the fluid. The evolution 



is similar to that of Penston (1969) and Foster & Cheva- 



lier (1993) for a collapsing pre-protostellar object, prior 
to what they call "core formation" (when a singularity 
forms) . 

For the run with a large number of initial Jeans masses 
(5.1Mj), significant turbulent motions are still present 



when the simulation is stopped. The fluid has condensed 
into a massive clump, surrounded by filamentary struc- 
tures, and many smaller clump-like objects. The individ- 
ual cores have a shallower density profile than those in the 
weak-gravity simulation, resembling that of a logatropic 
sphere (p oc r -0--^)) ( |McLaughlin fc Pudritz 1997| ). 
This suggests that they are still largely supported by 
turbulent motions. Similar density structures have been 



seen in regions of high-mass star formation (Van der Tak 
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et al. 1999; Van der Tak et ah 2000). In our simula- 



tions, the filaments and clumps appear to form simul- 
taneously out of the collision and merger of short shock 
arcs, rather than undergoing hierarchical collapse. How- 
ever, the shocks that do not merge do not contain enough 
mass to undergo gravitational collapse. 

We are currently performing and analyzing results 
from higher-resolution simulations. These will be presented 
in a future paper. 
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